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We study fermions in an electrically-probed and asymptotically anti-de Sitter Schwarzschild space- 
time which interact via novel chiral symmetry-preserving interactions. Computing the dual fermion 
two-point correlator, we show that these bulk interactions anisotropically gap Fermi surfaces of the 
boundary spectrum. Consequently, the interactions we devise provide holographic models for Fermi 
arcs seen ubiquitously in the pseudogap regime of the cuprates. Our interactions are modifica¬ 
tions of the chiral symmetry-breaking Pauli coupling, which has previously been proposed as the 
holographic realization of Mott physics. The onset of Mott insulation and pseudogap physics are 
respectively discussed in the context of bulk chiral and boundary parity symmetry breaking, and 
the Mott transition is interpreted as a deconfinement transition of non-Fermi liquid excitations. 


Superconductivity in the copper-oxide ceramics re¬ 
mains unresolved largely because of the unconventional 
electronic properties of the normal state. For example, 
when holes are doped into the copper-oxide planes, the 
metallic state that ensues is not characterized by a con¬ 
tinuous closed surface in momentum space as dictated by 
Landau’s paradigmatic theory of metals. Rather, the sur¬ 
face is truncated, forming what are referred to as Fermi 
arcs [1-5]. This stark deviation from the standard the¬ 
ory of metals can be viewed in one of two ways: either 
some type of order [6, 7] gives rise to a Fermi pocket 
with momentum-dependent spectral intensity that is van¬ 
ishingly small for some range of momenta, or the prob¬ 
lem is inherently rooted in strong coupling physics in 
which zeros of the single-particle electronic Green func¬ 
tion, caused by a divergent self-energy [8-12], are at the 
base of the vanishing spectral weight. These scenarios 
are distinguished based on their adherence to the Lut- 
tinger sum rule [13]. Within the former, quasiparticles 
exist but carry spectral weight too small to be measured 
experimentally on the ‘back half’ of the Fermi pocket, 
and so Luttinger’s rule is satisfied. In the latter, how¬ 
ever, whenever det Re G(a; = 0, k) = 0 the Luttinger rule 
is inapplicable [13-15]. Experimentally, the measured 
Fermi surface areas [2, 3, 5] at zero magnetic field violate 
Luttinger’s rule, a state of affairs which persists even at 
large fields [16]. The central problem of Fermi arcs thus 
appears to be elucidating how strong interactions per¬ 
sist from the Mott insulating state and partially gap the 
Fermi surface of the doped state. 

An ideal resolution of this problem would utilize a 
non-perturbative method to account for the strong in¬ 
teractions in the Mott state while delineating a mech¬ 
anism for Fermi arc formation. While Fermi arcs have 
been obtained numerically [12] and phenomenologically 
[10], they have evaded analytical methods that are valid 
in the strongly coupled regime. To address this short¬ 
coming, we utilize the gauge/gravity duality [17-19] (or 
‘holography’) — a method relating the physics of strongly 
interacting quantum systems to that of weakly interact¬ 


ing gravitational systems in one higher dimension — to 
investigate Fermi arc formation from non-Fermi liquid 
states. Prior holographic studies have generated Fermi 
arcs by anisotropically condensing fermions into p-wave 
[20] or d-wave [21] superconducting states, but this can¬ 
not describe cuprate physics because the cuprate arcs 
form at temperatures above the superconducting transi¬ 
tion. In contrast, we present a mechanism for arc forma¬ 
tion which does not utilize superconductivity, and there¬ 
fore represents a distinctly different state. The key result 
of this work is quite clear. States in which a gap forms 
without manifest symmetry breaking, hereafter referred 
to as Mott states, are realized at the boundary of asymp¬ 
totically anti-de Sitter geometries by bulk fermions un¬ 
dergoing chiral symmetry-breaking interactions. In con¬ 
trast, Fermi arcs obtain from the breaking of a discrete 
symmetry at the boundary, induced by chiral symmetry¬ 
preserving interactions in the bulk. We provide con¬ 
crete illustrations of the latter in bottom-up models with 
(2 -|- l)-dimensional and parity-broken boundary duals, 
and argue that boundary theories of chiral symmetry- 
invariant holographic fermions are most naturally inter¬ 
preted as two-fluid models undergoing momentum space 
confinement/deconfinement transitions. 

While there are several ways [19, 22-25] to implement 
the holographic program for fermion matter at finite den¬ 
sity, we pursue a bottom-up construction in which the ac¬ 
tion for a bulk (d -|- l)-dimensional gravitational system 
is supplemented with fermionic helds that source opera¬ 
tors at the d-dimensional boundary. Since this procedure 
provides only correlation functions for the boundary the¬ 
ory, there is considerable lee-way in choosing the fermion 
interactions in the bulk. Schematically our action 

‘5'bulk ‘^grav F F (1) 

will consist of gauge and gravity sectors, 5'gauge and S'grav 
respectively, with S'fer describing probe fermionic helds 
which source fermion operators in the boundary confor¬ 
mal theory. The only restriction on the gauge and grav¬ 
itational parts of the action is that they provide geome- 
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tries which asymptote to anti-de Sitter (AdS) spacetime 
at the boundary. In fact, one of the key conclusions of our 
work is that our results are independent of the detailed 
gravitational structure of the bulk, provided it contains 
electromagnetism and some form of a black hole. 

Our focus at the outset is the fermionic part of the 
action. Fermi to non-Fermi liquid behavior has pre¬ 
viously been shown [23, 24] to emerge from the sim¬ 
ple choice of the Dirac action minimally coupled to 
four-dimensional Reissner-Nordstrom-AdS (RN-AdS 4 ), 
in which bulk fermions "0 of mass m source boundary 
fermion operators of scaling dimension d/2 ± mL, with 
L the curvature radius of the asymptotic AdS geome¬ 
try. The two possible scalings of the boundary operators 
follow from independent prescriptions for identifying the 
operators’ sources and responses, referred to as standard 
and alternative quantization (resp. ±). This construction 
also produces gapped spectra as m is increased, but this 
has not been tied to Mott physics and the model’s sym¬ 
metry forbids descriptions of pseudogaps. Nonetheless, 
there are a number of non-minimal gauge interactions 
that can be added to extend this model, the simplest of 
which is the Pauli coupling, 

Sfer = [ i-ljj (0 — m — ipf) ip, (2) 


with p controlling the strength of a dipole interaction 
between the fermionic and Maxwell fields. Structurally, 
the dipole interaction provides shifts in fermion momenta 
that depend on the the boundary chemical potential; this 
charge scale shifts only the fermion frequencies in the 
minimal model. Tuning p from large negative (< — 1) 
to large positive (> 1) values (in the conventions of Ref. 
[26]) results in diverse phenomenology of the boundary 
fermions: the dominant low-frequency pole in their spec¬ 
trum passes through regimes of Fermi liquid-, marginal 
Fermi liquid-, and non-Fermi liquid-like scaling before 
reaching a gapped phase. While the gapped structure 
was thought [26, 27] to obtain from a vanishing of the 
quasiparticle residue, in actuality it arises from an ex¬ 
act pole-zero duality within the diagonal blocks of the 
boundary fermion propagator. 


Gii{uj,k;m,p) 


1 

Gii{u},-k;-m,-p)’ 


(3) 


of the spectral weight is due to zeros and the bulk Pauli 
coupling in RN-AdS 4 and electrically-probed SS-AdS 4 
mimics Mott physics. 

There is a subtlety, however, in the Mott gap formed 
from the Pauli interaction. Although the symmetries of 
the boundary spectrum are preserved, the Pauli term has 
a non-zero anticommutator with the generator of chiral 
rotations, {r“*’,r 5 } ^ 0. That is, chiral symmetry is 
broken in the bulk. This is not entirely surprising since 
chiral symmetry breaking is a typical mechanism for the 
generation of mass. However, more than the static break¬ 
ing of this symmetry is relevant here. What is crucial to 
note is that the Pauli term only generates a gap for suf¬ 
ficiently large and positive values of p; there would be no 
such restriction if the gap were attributable merely to loss 
of chiral symmetry. Because the Pauli term changes the 
scaling dimension of the dual boundary operators and in¬ 
creasing the exponent converts pole singularities to zeros, 
it is ultimately the dynamical breaking of chiral symme¬ 
try that is the root cause of the gap. The same reasoning 
applies to the gap generated by the Dirac mass m in the 
minimal model, though in that case gapping ensues with¬ 
out alteration of the coupling between bulk fermions and 
the boundary charge scale. Though the dimensionality of 
the bulk and boundary differ, they share time coordinates 
and fermion charges and so must share time reversal and 
charge conjugation symmetries. Bulk chiral symmetry 
breaking must then be reflected in discrete symmetries of 
the boundary. Extrapolating to flat-space lattice models, 
where chiral symmetry is a combination of particle-hole 
and time reversal symmetry, we may take the holographic 
results to mean that the generation of Mott-type gaps re¬ 
quires the breaking of one of these symmetries. 

Because the pseudogap is not a completely gapped 
phase, we expect holographic models for Fermi arcs to 
preserve bulk chiral symmetry. Exploiting the pole-zero 
duality inherent in the Pauli construction, we anticipate 
that a bulk interaction which couples fermion momenta 
to different signs of the Pauli term should generate both 
poles and zeros and hence Fermi arcs. Consequently, we 
propose the following bulk fermion action to model Fermi 
arcs, 

<5'fer = J dP^Xy/^ iip (^0 - m - ipiT0 + p 2 (n ■ ip. 

(4) 


first shown for RN-AdS 4 [28] and later for electrically- 
probed Schwarzschild-AdS 4 (SS-AdS) [29]. The in¬ 
verse relationship [23] between Gn(uj,k\m,p) and 
Gii{uj, —k] —m, —p), ultimately rooted [29] in the two 
quantizations for holographic fermions, was not exploited 
until it was realized [28] that the boundary spectrum 
must solely exhibit zeros for large positive p, as it con¬ 
tains only poles at large negative p. It is well-known 
[8, 9, 11-13, 30] that the Mott gap requires zeros of the 
single-particle propagator. Consequently, the vanishing 


The non-minimal interactions tuned by pi and p 2 dif¬ 
fer from the dipole interaction of (2) through the pres¬ 
ence of the matrices F = rM-(h • F) and n ■ F, which 
restore bulk chiral symmetry while breaking rotational 
and Lorentz symmetries of the boundary theory. The loss 
of boundary rotational invariance is necessary to model 
anisotropic phenomena such as Fermi arcs and is charac¬ 
terized in our model by the breaking of parity along the 
unit vector n.[37] For simplicity we examine the interac¬ 
tions individually. 



3 


Beginning with the pi interaction, we may naively infer 
the effects of F from our understanding of the pole-zero 
duality and the Pauli coupling model. From the work 
of Refs. [26, 29] we know that, with m = 0, the Pauli 
model produces sharp Fermi surfaces at large negative 
p and gapped spectra at large positive p. In terms of 
the two diagonal entries of the boundary fermion Green 
function k), these Fermi surfaces manifest in two 

poles: one of Gn at (w = 0, |kj?|), and one of G 22 at 
(w = 0, —|ki?|), with the Fermi momentum k^? depend¬ 
ing on the background geometry. The pi interaction 
enters the bulk Dirac equation in the same way as the 
Pauli coupling along the h-momentum axis — preserv¬ 
ing the general pole/zero structure of Gn and G 22 in 
this frame — but inverts the sign of p in the Gn block. 
Consequently, in the n frame the pole that was once 
present kp becomes a zero at —kp, coincident with 
the pole still present in G 22 , and when G 22 ( 0 , kp\p) and 
Gii(0 ,/cf; —p) have negligible spectral weight the Fermi 
surface gaps dX kp. Conversely, the p 2 interaction is a 
less straightforward modification of the dipole coupling, 
and therefore does not manifestly engineer pole/zero co¬ 
existence in the boundary spectrum. [38] Its virtues are 
preservation of bulk chiral symmetry and omission of the 
radial boost generator F-F- present in the pi interaction; 
the former property is demanded by our line of reasoning, 
while the latter gives the interaction a form more aes¬ 
thetically natural from the perspective of the boundary, 
where radial boosts lack straightforward interpretation. 

While the above considerations provide some intuition 
about the spectrum in the h frame, the model’s lack 
of rotational symmetry makes it difficult to infer spec¬ 
tral properties at general momenta. To remedy this we 
have numerically computed the boundary fermion spec¬ 
tral function in planar SS-AdSi with a Maxwell probe. 
Qualitatively identical results are found in extremal RN- 
AdS 4 . Unnormalized spectral densities of the boundary 
fermions at the small frequency ui = 10“^ -I- iS, for four 
values of the pi interaction and two of p 2 , both with 
n = X, are presented in Fig. 1. A small broadening 
factor 6 = 10“® was introduced to resolve poles in the 
retarded Green function. More information about the 
geometry and computation of Green functions, including 
parameters used to generate Fig. 1, may be found in the 
appendix. 

In the absence of the non-minimal interactions, 
electrically-probed SS-AdS 4 hosts a highly-broadened 
Fermi surface in the boundary dual. Once either inter¬ 
action is switched on, the previously-discussed gapping 
process ensues: spectral weight on the right half of the 
Fermi surface is suppressed, while that on the left half 
is enhanced. This trend continues as pi and p 2 are in¬ 
creased, leading to an unmistakably arc-like spectrum as 
Fig. 1 reveals. If either coupling constant is instead tuned 
to negative values, the spectral suppression and enhance¬ 
ment occur on the left and right halves of the Fermi sur¬ 


face (respectively). Thus, our bulk interactions neces¬ 
sarily provide gapless spectra in the boundary dual, in 
contrast with the chiral symmetry-breaking interactions 
that have previously been studied. Note that the cen¬ 
ter of mass momentum produced by non-zero pi, which 
is significantly less drastic for non-zero p 2 , may be can¬ 
celed by introducing a second, independent flavor of bulk 
fermions that experience the same interaction with cou¬ 
pling —pi. In addition, our mechanism generates arcs 
without the aid of superconductivity [20, 21], and hence 
could provide a framework for understanding the emer¬ 
gence of arcs in the cuprates. 

From the bulk perspective the arc formation is most 
easily understood in terms of bulk fermion orbits. Tech¬ 
nically, poles and zeros of boundary fermion propaga¬ 
tors reflect the linear independence of fermion sources 
and responses; these objects are identified with the 
boundary values of the two (two-component) F- eigen- 
spinors, and the propagator is defined as the transfor¬ 
mation matrix that relates them (see appendix). Yet 
the bulk is nothing but an on-shell scattering problem 
for Dirac particles in a radial potential encoding grav¬ 
ity and electromagnetism. [39] Numerical solutions of the 
Dirac equation (with in-falling boundary conditions) re¬ 
veal that the F- eigenspinors assume non-normalizable 
scattering states for most boundary wave modes. At 
certain frequencies and momenta, however, one of the 
spinors (or more generally one of its components) can 
enter a normalizable bound state; this is precisely where 
poles and zeros appear in the boundary spectral func¬ 
tion. Holographic Fermi surfaces are attributed [23] to 
the ‘source’ spinor occupying a bound state. Similarly, 
zeros appear at frequencies and momenta where the ‘re¬ 
sponse’ spinor is bound to the black hole. The rela¬ 
tionship between the pole-zero duality and the choice 
of quantization scheme naturally follows. Our model (4) 
thus achieves a coexistence of poles and zeros, and ex¬ 
hibits Fermi arcs, by mediating bound states for both 
types of fermions via its non-minimal interactions. 

Understanding the arc from the boundary perspec¬ 
tive is a more difficult task, as the field theory mecha¬ 
nisms behind holographic Fermi surfaces are under de¬ 
bate (competing interpretations are discussed in Refs. 
[31] and [32]), while those behind holographic zeros have 
not been investigated. We can make some interpretive 
progress, however, by appealing to holographic general¬ 
izations of Luttinger’s rule [33]. These generalizations il¬ 
lustrate how on-shell Dirac fermions — i.e., bulk charges 
external to a black hole — provide a traditional Luttinger 
count by adding the volumes contained by Fermi surfaces 
to the boundary charge. In contrast, charges contained 
by a black hole provide an explicit deviation from the tra¬ 
ditional count. But there is a critical oversight in these 
works as they ignore the possibility of zeros. The Pauli 
model provides many instances of zero surfaces in the 
boundary spectrum, and by the arguments here zeros 
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FIG. 1. Boundary fermion spectral functions A{uj, k^, ky), at a; = 10“^+ 15, 5 — 10“®, in the Fermi arc model (4) with couplings 
of pi = 0, 1, 2, 4 with p 2 = 0 (panels (a)-(d)), and with p 2 = 0.2, 0.4 with pi = 0 (panels (e)-(f)). Other bulk parameter 
choices are listed in the appendix. 
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are attributable to charges outside the black hole. Thus, 
even in the holographic context, zeros represent charged 
degrees of freedom which lack structure in momentum 
space. Our Fermi arc model then implies that the pseudo¬ 
gap is a two-fluid state composed of momentum-confined 
and momentum-deconfined charges. Building a field the¬ 
oretic understanding of the latter would not only be in¬ 
teresting in its own right, but may also elucidate how 
pseudogaps intervene the development of non-Fermi liq¬ 
uids from Mott insulating states. 

Given the central role played by gravity and electro¬ 
magnetism in the bulk, it is natural to wonder how back- 
reaction of fermions will affect results from probe mod¬ 
els. The key questions concern how bulk geometries and 
fermion orbits are altered by backreaction. Fermion or¬ 
bits are observed widely enough in bottom-up and top- 
down models for us to expect their presence in generic 
AdS-black hole geometries, but it is not obvious that 
backreacted systems will support the a; = 0 bound states 
necessary for boundary zeros and Fermi surfaces. What 
is encouraging, however, is that recent work on the ex¬ 
istence of Fermi surfaces in top-down constructions [34] 
can be understood entirely in terms of a competition be¬ 
tween gravity and electromagnetism. In the cases studied 
[34], Fermi surfaces are always observed when there are 
bulk fermions with (positive) charge appreciably greater 
than their mass, provided the black hole hosts the corre¬ 
sponding electric field. In marginal cases where fermion 
charges and masses are comparable, the existence of a 
Fermi surface is contingent upon the presence of a posi¬ 
tive Pauli interaction. Hence, investigating backreacted 
geometries should reveal whether or not the presence of 
zeros leads to an increase in the charge behind the hori¬ 
zon as a result of infall of fermion bound states. Such 
infall would be consistent with Mott insulation arising 
from deconfined charges behind the horizon. 
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Appendix: Model details and calculation of 
correlators 

We study bulk fermions with non-minimal interactions, 
described by the Lagrangian 

Aer = -m- ipiVf + p 2 {n ■ Ip, (5) 


with F = rM^(n • f) and f = (F^^, Fit). We take n = x 
henceforth. The covariant derivative and Maxwell tensor 
may be written 


f -f'MAT, 


( 6 ) 


with the (inverse) vielbein, uj'j^ the spin connection, 
and = i[r“,r^]. Our index conventions use capital 
Roman letters for bulk coordinates M, N ■ ■ ■ = {t,x’‘,r} 
and lower-case Roman letters for tangent space coordi¬ 
nates a,b - ■ ■ = {t,x^,r}. 

For the geometric background we choose 
Schwarzschild-AdS in (d -|- 1) = 4 dimensions. 

Parametrizing in the Poincare patch, the line ele¬ 
ment may be written 

ds2 = ^(-/(r)dt2 + dx2) + ^^. (7) 


We take the black hole to have unit (dimensionless) mass, 
so the emblackening factor is given by 

/(r) = l-(^)', ( 8 ) 

with ro denoting the horizon radius. The temperature of 
the boundary theory is then T = 3ro/ 47 rL^. We give the 
Maxwell probe the form familiar from RN-AdS 4 , 

A = dt, Ai = Qro/L^, (9) 

and treat Q as a tuning parameter for the boundary 
chemical potential pL. When we refer to electrically- 
probed SS-AdS in the main text, we mean the metric 
(7) with the probe Maxwell field (9). The plots in Fig. 
1 were generated with Q = -s/S, rg = L = q = 1, and 
m = 0. 

To evaluate the Dirac operator it is convenient to 
Fourier transform the bulk spinor in the boundary co¬ 
ordinates and scale it by the factor to elim¬ 

inate the spin connection. Hence we take ip(r,x) ^ 
' 0 (r, with k = (w, k). Ghoosing the fol¬ 

lowing basis for the Dirac matrices, 


rr 

pi 


0 \ pt = 0 \ 

'v 0 -as) ’ 0 zaiyl ’ 

f-a2 0 \ pz ^ 0 cr2\ 


( 10 ) 


and rescaling the non-minimal couplings as pi —)■ ptL, 
the Dirac equation for the scaled fields yields two coupled 
equations. 
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(ui + qfi li¬ 


ra 


- {-lycfi -ki 

\ r 


)) V'i - 0-3 + (-l)-^P2Ai^) 1p3 

^ yij + crifc2V’i- 


( 11 ) 


Here we have expanded the spinor as ijj = (' 0 i,V' 2 )^, 
with the two-component spinors ijjj reflecting the block 
structure of (10). The final factor of ipi should strictly 
take i y j. For a more detailed derivation, see Ref. [26]. 

The retarded Green functions of boundary fermion 
operators are realized by asymptotic solutions to 
the bulk Dirac equation, subject to in-falling condi¬ 
tions on the bulk spinors at the black hole horizon. 
Asymptotically the solutions to (11) are '0j(r, fc) = 
{bj{k)r~'^^,aj{k)r'^^Y', with j G {1,2}. In the mass 
window |toL| < 1/2 we are free to choose the aj or the bj 
as the sources for fermion operators in the dual theory. 
Though there is no distinction between these “quanti¬ 
zation” schemes in the zero mass case of the minimal 
model, we make the standard choice of A = ( 01 , 02 )”^ 
as the source and B = ( 61 , 62 )^ as the vacuum expecta¬ 
tion value of a boundary fermion operator. Generally the 
source and vev are related through a linear transforma¬ 
tion, B = iSA, from which the Green may be computed 
via [35] 

G(w,k) = -^57^ ( 12 ) 

with 7 - a boundary Dirac matrix ( 7 - = iai in our basis). 
The causal structure of this Green function is determined 
by boundary conditions on the bulk spinors at the black 
hole horizon, with in-falling conditions providing the data 
for retarded correlators. 

Though it is possible in principle to compute the trans¬ 
formation matrix S through numerical integration of the 
Dirac equation, in practice it is often easier to work with 
a set of first-order, non-linear evolution equations for the 
Green function components. To implement this proce¬ 
dure, we first expand the Dirac equation about the hori¬ 
zon to determine a basis of in-falling states for the bulk 
spinors. The near-horizon expansion gives 

where Cj = luL^/ r^d. Writing the two-component bulk 
spinors as in-falling solutions are 

given by the eigenvectors of diag(z(T 2 , * 0 - 2 ) with eigenvalue 
—i , or 

= (i, 1,0,0)^ , /3l=i, a\ = 1, j3\ = a\ = 0, 

(14a) 

= (0,0, i, 1 )^ , /?" = = 0, jdiy = i, = 1. 

(14b) 


The components aj and j3j act respectively as sources 
and responses for fermion operators at the conformal 
boundary, with the Green function given by (12). The 
data from both sets of boundary conditions can be en¬ 
coded in a single equation by defining the matrices [36] 




Y = GZ, (15) 


such that G asymptotically realizes the boundary Green 
function (12). Taking a derivative of the third equality 
and making use of the Dirac equation (11), we derive the 
following evolution equation for G, 

-^JJ{r)drG = M+ + GM + MG-GM_G, (16) 
where the matrices M± and M are 

±v±{r) - koo ky 

ky ±v^ (r) -I- ka: 

M=-( ^ 

r \ 0 -p 2 /iroy'’ 

and 

^^±(r)= ^^(a;-fgAt(r))±pi/i—. (18) 

The in-falling conditions (14) provide the initial condition 




lim G = 

r^ro 



(19) 


which allows numerical integration of (16). 
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